Fate of non-Fermi liquid behavior in QED 3 at finite chemical potential 
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The damping rate of two-dimensional massless Dirac fermions exhibit non- Fermi liquid behavior, 
oc e 1//2 , due to gauge field at zero temperature and zero chemical potential. We study the fate 
of this behavior at finite chemical potential. We fist calculate explicitly the temporal and spatial 
components of vacuum polarization functions. The analytical expressions imply that the temporal 
component of gauge field develops a static screening length at finite chemical potential while the 
transverse component remains long-ranged owing to gauge invariance. We then calculate the fermion 
damping rate and show that the temporal gauge field leads to normal Fermi liquid behavior but 
the transverse gauge field leads to non-Fermi liquid behavior oc e 2//3 at zero temperature. This 
energy-dependence is more regular than oc e 1 / 2 and does not change as chemical potential varies. 
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I. INTRODUCTION 

The damping rate of fermions due to interaction with 
gauge fields is a physical quantity of broad interests. 
Studying this quantity can help us to judge whether an 
interacting fermion system displays non-Fermi liquid be- 
havior or not. According to Landau, for any normal 
Fermi liquid to be stable, the fermion excitations must 
have a sufficiently long lifetime, which means that the 
fermion damping rate should vanish faster than energy 
e does as e — ¥ 0. At the low-energy regime, the fermion 
damping rate can be written in the form ImS(e) oc e z . 
The system with exponent z > 1 is a normal Fermi 
liquid, while system with z < 1 corresponds to a non- 
Fermi liquid. In conventional metals, the Coulomb in- 
teraction between electrons is always statically screened 
and can only lead to normal Fermi liquid behavior. Since 
the work of Holstein and coworkers^, it has been known 
that the unscreened gauge field can give rise to non- 
Fermi liquid behavior. The unusual, non-Fermi liquid 
like, fermion damping rate caused by gauge field has at- 
tracted great attention in the past twenty years because 
an emergent gauge field is found to play important roles 
in a number of strongly correlated electron systems^—. 
In particle physics, this problem has also been discussed 
extensively in various gauge theories, including four- 
dimensional QCD 9 and four-dimensional QE D 10 ' 11 . 

Here we are particularly interested in the unusual 
properties of massless Dirac fermions. In some planar 
correlated electron systems, including d-wave high tem- 
perature superconductor— and graphenei^, the valence 
band and conduction band touch only at discrete Dirac 
points. This is illustrated in Fig. 1(a). The states in 
the lower valence band are fully occupied, while those 
in the upper conduction band are fully empty. The 
low-energy fermions excited from the lower band have 
a linear spectrum and hence can be described by mass- 
less Dirac fermions, which satisfy the relativistic Dirac 
equation&i£. The interaction of massless Dirac fermions 
with an abelian gauge field in two spatial dimensions 
defines the three-dimensional quantum electrodynamics 



(QED3). This field theory is usually studied by parti- 
cle physicists as a toy model of QCD4 since it is known 
to exhibit dynamical chiral symmetry breaking— and 
confinement^. In condensed matter physics, with proper 
modification, it serves as an effective low-energy theory 
of high temperature superconductors^— and some spin- 
1/2 Kagome spin liquids^. In the realistic applications, 
chiral symmetry breaking and confinement in QED3 cor- 
respond to the long-range Neel order in two-dimensional 
quantum Heisenberg antiferromagneti^. 

Recently, we studied the damping rate of massless 
Dirac fermions due to gauge field in QED 3 in the ab- 
sence of chiral symmetry breaking and confinement 20 . 
It diverges at both zero and finite temperatures when 
it is calculated using the straightforward perturbativc 
expansion approach. Once the fermion damping effect 
and the dynamical screening effect of gauge field are 
self-consistcntly coupled, the fermion damping rate is 
then well-defined and behaves as ImE(e) oc e 1 / 2 at zero 
temperature^. This damping rate vanishes slower than 
e does in the e —¥ limit and thus displays non-Fermi 
liquid behavior. 




FIG. 1: (a) The half-filling state with lower band being fully 
occupied and upper band fully empty, (b) At finite chemical 
potential fi, the Fermi surface is finite. 

This result was obtained in the so-called half-filling 
ground state depicted in Fig. 1(a). In realistic condensed 
matter systems, the fermion density already can be con- 
tinuously turned, either by chemical doping 8 or by ad- 
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justing a bias voltage^. When the fermion density grows 
from the Dirac point, a small but finite Fermi surface will 
emerge, as shown in Fig. 1(b). To describe this process, 
the commonly used strategy is to introduce a chemical 
potential fi, which defines the energy difference between 
the new Fermi surface and the Dirac point. The systems 
at zero and finite /i may have quite different properties. 
Indeed, there might be a quantum phase transition as 
the chemical potential varies 2 ^, with [i = being the 
quantum critical point. 

When the fermion density becomes sufficiently large, 
the interacting fermion system finally develops a large 
Fermi surface with low-energy fermionic excitations sat- 
isfying the non-relativistic Schrodinger equation. His- 
torically, before the energy gap of high temperature 
cuprate superconductor was confirmed to have a d-wave 
symmetry, the fermion-gauge system with large Fermi 
surface had been studied intensively 2 ^ - — . The damp- 
ing rate of non-relativistic fermions was calculated by 
various methods, including straightforward perturba- 
tion expansion^, renormalization group approach^, and 
Eliashberg theory 7 . Most of these studies found that 
ImE(e) oc e 2 / 3 at zero temperature. The exponent in the 
energy dependence of damping rate is quite different from 
that in QED3 of massless Dirac fermions at \i = 0. The 
difference is presumably owing to the difference between 
a large Fermi surface and discrete Fermi points. 

As fi grows from fi — to a large value, the fermion 
damping rate will undergo a crossover from cx e 1 / 2 to 
oc e 2 / 3 . A question naturally arises: Does the exponent 
z appearing in e z vary continuously from 1/2 to 2/3 or 
change abruptly at some critical value \i c l 

In this paper, we study how the fermion damping rate 
varies with growing chemical potential fi by calculating 
the /^-dependence of fermion self-energy. As usual, the 
gauge field is decoupled to longitudinal and transverse 
components. From the vacuum polarization functions at 
finite /z, we know that the longitudinal component be- 
comes short-ranged (massive) due to static Debye screen- 
ing effect, but the transverse part remains long-ranged 
(massless) because of the gauge invariance. In this case, 
the transverse component of gauge field dominates and 
should be able to produce non-Fermi liquid behaviors. 
However, with increasing fermion density, the dynamical 
screening effect becomes stronger and hence may lead to 
less singular behavior than that at small ji. 

After explicit computation, we found that the trans- 
verse damping rate of Dirac fermions at zero temperature 
behaves as ImSx(Ck) oc A 1 " 1 ^ 3 ^ 3 ' while the longitudinal 
contribution is ImSL(Ck) oc (££//•*) m (£k/M)> where £k is 
the fermion energy in the on-shell approximation. Thus 

the total fermion damping rate is oc M ^ 3 £k^ 3 j which 
is certainly non-Fermi liquid behavior. We also consid- 
ered the fixed momentum approximation and obtained 
the same results, i.e., the fermion damping rate behaves 
as oc //~ 1 / 3 e 2 / 3 at zero temperature. 

These results imply that the fermion damping rate sud- 
denly becomes oc e 2 / 3 from oc e 1 / 2 once the chemical 



potential \i departs from zero. As /x grows, the energy- 
dependence of fermion damping rate does not change 
but its coefficient decreases. Therefore, although the 
Dirac fermions are always not well-defined in the sense of 
Landau quasiparticle, their lifetime increases slowly with 
growing fermion density. 

The paper is organized as follows. The Lagrangian and 
some relevant quantities are defined In Sec. II. The full 
expressions of polarization functions from massless Dirac 
fermions at finite chemical potential are calculated in Sec. 
Ill and the fermion damping rate is calculated in Sec. IV. 
We summarize the results and briefly discuss the physical 
implications in Sec.V. The polarization functions at zero 
temperature are shown in the Appendix. 



II. LAGRANGIAN AND FEYNMAN RULES OF 
QED 3 AT FINITE CHEMICAL POTENTIAL 

We start from the following general Lagrangian of 
QED 3 

N 



1=1 



/i — ieao — ia ■ (d — iea)) — —F' ' . 

(1) 

This is the general Lagrangian for QED defined in (2+1)- 
dimensional space-time. As a well-defined relativistic 
quantum field theory, there is surely an explicit Maxwell 
term for the gauge field. As discussed in the Introduc- 
tion, this is a very interesting and widely studied field 
theory in the context of particle physics. When applied 
to strongly correlated electron systems, it usually needs 
to be modified properly. If the effective QED3 theory is 
derived by considering the phase fluctuations in under- 
doped high temperature superconductors, then normally 
the Maxwell term is present. However, if the effective 
QED3 theory is obtained by the slave-particle treatment 
of t-J model, there is no Maxwell term in the Lagrangian 
and the gauge field can have its dynamics only after inte- 
grating out the matter field o 15 ' 16 . We will first consider 
the general action with the Maxwell term and briefly dis- 
cuss the case without such term at the end of Sec. IV. 

Here, we adopt the two-component representation of 
spinor field, with a being the Pauli matrices. The Dirac 
fermion flavor N is taken to be large so that we can use 
the 1/N expansion. The theory is defined at finite chem- 
ical potential /j,. The aim of this paper is to study how 
the fermion damping rate depends on [i. For simplicity, 
we take h = c = ks = 1 throughout the whole paper. 

At finite temperature, the Matsubara propagator of 
massless Dirac fermion is 



G (ie n ,k) = - 



1 



(2) 



is n + fi — a ■ k 

where e n = (2n + 1)ttT with n being integer. After ana- 
lytic continuation, the retarded propagator is 

1 



G (e,k) 



fx — a • k + id 



(3) 
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At finite temperature, the temporal and spatial compo- 
nents of gauge field decouple and now it is convenient to 
work in the Coulomb gauge fcjdi = 0. In the imaginary 
time formalism, the propagator for the gauge field can 
now be written as 



Aw («w m ,q) 



Dij (iUmil) 



1 



|q| 2 + n 00 («w m ,q) ' 



(4) 



1 



n_L (iu) m ,q) 



(5) 



where ui m — 2mnT for bosonic modes with m being in- 
tegers. The vacuum polarization functions n o(w m ,q) 
and Hj_(uj mi q) come from the one- loop bubble diagram 
of Dirac fermions to the leading order of 1 /N expansion. 
In particular, the polarization function appearing in the 
spatial component is given by 



uj- 



n_i_(«w m ,q) = Uu(iu) m ,q) Sn 00 (iu; TO ,q) 



(6) 



n 00 (iu; m ,q) 



ILii(iu m ,q) 



The functions IIoo(iw m , q) and ILa(iu) m , q) are defined as 

9 m ^-^ f d 2 k 
—Ne T} / - — ^ 

xTr[G (je„,k)Go(i£„ + iw m ,q + k)], (7) 

x Tr [o-i G (i£„, k.)aiG Q (ie n + iu> m , q + k)]. 

(8) 

When we calculate the fermion damping rate, we need 
the real and imaginary parts of the retarded polariza- 
tion functions. They can be obtained by straightforward 
computation, which will be given in the next section. 

The fermion damping rate can be calculated by the 
standard finite temperature field theory technique^ 2 -. To 
the lowest order of 1/N expansion, the self-energy of 
Dirac fermion is given by 



where 
S L (i£: n ,k) 



E(ie n , k) = E L (ie„, k) + E T (ie n , k), 

d 2 q 



(9) 



1. 



2 T]T -Tr 



{2-Kf 



xG (is n + iw m ,k + q)£)oo(iw m ,q) , (10) 



E T (fe n ,k) = e 2 Tj2^ 



_cfq_ 

(2tt)2 



XCTiG (ie„ + iuj m , k + q)a j D lj (iuj m , q) 

(11) 

are the contributions from the longitudinal and trans- 
verse components of the gauge field, respectively. The 



damping rate of massless Dirac fermion will be obtained 
by making analytic continuation, ie n — > e + iS, as 



E(e,k) = E L ( £ ,k) + E T ( e ,k), 



(12) 



and then taking the imaginary part, ImE(e,k). 

When the Fermi level lies exactly at the Dirac point, 
the states below the point are all occupied while the 
states beyond it are all empty (see Fig. 1(a)). In this state, 
the chemical potential is usually defined as zero, fj, = 0. 
In a previous paper, we studied the Dirac fermion damp- 
ing rate and found that it behaves as ImE(e) cx e 1 / 2 at 
zero temperature, which is a typical non-Fermi liquid be- 
havior. Once the fermion density increases starting from 
the Dirac point, the system develops a finite chemical 
potential \i (Fig. 1(b)). Now the system has a finite but 
small Fermi surface. The density of states of fermions at 
the Fermi level has a finite quantity. Therefore, at finite 
/j,, the gauge field may lead to very different behaviors 
for the fermion damping rate. 

In order to know how fermion damping rate varies with 
fi, we will explicitly calculate the fermion self-energy. 
To this end, we first calculate the polarization functions 
n o(iw m ,q) and Uu(iu m ,q). 



III. COMPUTATION OF POLARIZATION 
FUNCTIONS 

The polarization functions contributed by the massless 
Dirac fermions deserve careful exploration since they de- 
termine or are directly related to many important phys- 
ical quantities. For instance, the dynamical screening 
effect of collective particle-hole excitations on the gauge 
or Coulomb interaction between Dirac fermions can only 
be studied by the polarization functions. Physically, such 
effect describes the damping of gauge boson in the many- 
body background composed of massless Dirac fermions. 
In addition, according to the Kubo formula in transport 
theory, various conductivities are all given by their corre- 
sponding current-current correlation functions, which in 
form are analogous to the polarization functions. 

In this section, we briefly outline the computational 
steps and present the complete expressions for polariza- 
tion functions IIoo and LT_l in the presence of finite chem- 
ical potential at both zero and finite temperature. The 
analytical expressions will be useful to any work that 
relies on the properties of polarization functions of two- 
dimensional Dirac fermions. 



A. Temporal component IIoo(w,q, T) 

To calculate the temporal component of polarization 
function Hoo(ioj m , q), we first introduce the spectral rep- 
resentation 



n c m f +00du 
G {ie n ,k) = — 

J-oo t 



duii Im [Gq(wi, k)] 



(13) 



and then sum over the frequency, which yields 



n o (iuj m , q) = -Ne 



d 2 k 



rTr 



^ilm[G ( Wl ,k)] / ^Im[G ( W2 ,k + q)] 



(27T)2 

n F (u>i) - n F (LL> 2 ) 

UJl — U>2 + iuj m 

It is convenient to make the analytic continuation ioj m — > uj + iS at this stage: 

1 1 „ 1 



UJl — UJ2 + i^ri — UJ-2 + W + iS UJl — UJ2 + UJ 

The imaginary part of the retarded polarization function is 

j2i, r r+°° a..,. r+°° 



ittS^oji — u>2 + UJ )- 



— — Tr / -^Im[Go(wi,k)] / -^Im [Go (fa*, k + q) 
x [np(ui) — rip (fa*)] $ (fa>i — fa* + w) . 



Here the imaginary part of retarded fermion Green function is given by 
Im[G (w,k)] = Im ' 

= Im 



-7TS, 



uj + [i — a ■ k + iS 

uj + [i + a ■ k 
(uj + [if - |k| 2 + jsgn (uj + n) 6 
gn(fa) + fi) (uj + fi + a ■ k) 5 ((uj + fi) 2 - |k| 2 ) 



— 7rsgn(u; + fi) (uj + fi + a ■ k) —— ■ [S (ui + fi + |k|) + S (ui + fi — |k| 

Z K 



After tedious computation, we finally have 



1+e 2T 



Imlloo (w,q,T) = < 



when \uj\ > |q| 



l + e 2T 



-2o M q|.i + |^|-2e. M 

l + e 2T 



when a; I < |q| 



We now calculate the real part of temporal component of polarization function. The whole computation 
more complicated than the imaginary part. From Eq. (|14[) and Eq. (|15|) . we have 



Tr / ^llm[G ( Wl ,k)] / ^Im[G (fa*,k + q )] 
- n F (uj 2 ) 



UJl — UJ2 ~\~ UJ 

Ne 2 I ^Tr 



r d 2 k r /• +00 ^ 1 /• +00 du; 

/ — — -Tr / Im Go w h k % Wl P — 



efcj 2 Im [G (fa*,k + q)] 



-Ne 2 



' — oo 
c+oo 



f d 2 k T r 00 duj 2 ,„ , , S1 , . /" +00 dfa) 
/ — -Tr / — ±Im G fa*,k + q )}n F (cj 2 )P / — 



UJl +UJ-UJ2 

+°° dui Im[G (wi,k)] 



UJ2 — UJ — UJl 



Using the Krames-Kronig relation 
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the above expression can now be converted to 

/<i 2 k r+ ca 
-—2 / ^n F {u l )T^[lm[G Q {^ x ,V)]Re[G {u l +u,V + ^]] 



— OO 

/ ^^(wOTrpmlGoK^ + q^RelGoK-w^)] 



+Ne 



where 



(21) 



RcG (lo, k) = (oj + fx + a ■ k) P 



(o; + M ) 2 -|k| 2 ' 



From this equation, we obtain the following expression: 



Ne 2 T\n( l+e~ 



Ne 2 kl 2 r + °° dx\/x 2 - 1 



±1 "I 2tt 
1 



iVe 2 /"+ 00 
Ren 00 (w, q, T) = -— y rfjkj + <{ 



|q|r-|^||-2c. M q| J + |^|-2a M 

1+e 2T 1 + e 2T 

when \u\ > |q|, 
when |w| < |q|. 



"1 

5a,l 



1+e ST" 



Notice there appears a divergent term 



2 /-+0O 



-^Singular 



Ne 2 
~^ Jo 



d|k|. 



(22) 



(23) 



(24) 



To remove this divergence, here we employ the regularization scheme that was proposed in=,. This scheme states that 
the gauge field must remain massless so that it should satisfy 



IV (w -> 0, |q| -> 0,/i = 0,T = 0) = 0. 
Now the polarization function can be re-defined as 

IV(w,q,T)-IV (w->0,|q| ->0,/i = 0,T = 0). 
After this regularization, the singular term can be simply dropped from Relloo (w, q, T). 



(25) 
(26) 



B. Transverse component II x (o;,q, T) 



Proceeding as we have done in the above, we found that the imaginary and real parts of Ha (w,q, T) have the 
expressions: 



E a=±1 \ -sgn (w) ^ 7 J^ J\ dxVT^ 



q| 3 ;+|^|-2c. >J 

1+e 2T 



== q ' .r ■• ■ I -u' I — 2a- /./. 

1+e 2T 



Inalljj (w,q,T) = < 



when | a; | > |q| 



(27) 



q.' .-■ 2../, q|x+M-2aM 

1+e 2T 1+e 2T 



|q|a:-M~2Q M q | i + | u | - 2c M 

1+e 2T 1+e 2T 



when \u)\ < |q| 
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2c M q. r -|^|-2a)i 

1+e 2T 1+e 2T 



q .r ,c 2m /i | q | ,x -f- 1 oj | -2a(i 

1+e 2T 1+e 2T 



Ren« (w,q,T) = < 



when a; > |q|, 



V J jv£ Iql- f 1 w x yi^2 



1+e ?r 



|q|x + |^||-2a t i 

1 + e 2T 



when |w| < |q|. 



According to Eq.©, the retarded transverse polarization function is decomposed as 



(28) 



U ± (lj, q, T) = Iki{u, q, T) + ^n 00 (w, q, T), 



(29) 



which can be written more explicitly as 



ImII+(w,q,T) = lmU ll (u,q,T) + — Imn o(w,q,T), 



ReIl+(w,q,T) = Ren il (w,q,r) + ^Ren o(w,q,T). 



(30) 
(31) 



Using the results presented above, it is easy to get that 
Imn±(o;,q,T) = < 



S a ,1 



|q|i+|^|-2o.^ 

1+e 2T 



£ Q=±1 sgn(.)^^^2 j+oo dx _s_ 

1 



|q|s+|"|-2°H 

1+e 2T 



|q|a;- |oj I -2q^ 
.1 + e 2T 



when \uj\ > |q| 



when \u\ < |q| 



(32) 



Ne 2 T\n 



a=±l "I 2tt 
1 



(l+e^) 



|q| 2 8^ 
l 



||q|x- . |-2 Q)J q|;l+|^|-2c M 

1+e 2T 1+e 2T 



Rcn ± (a;,q,T) = < 



when \w\ > |q| 



x <5 Q ,i 



Afe 2 Tln 



(l+e^) 



2tt 



|q| 2 + ^ 



1+e 2T 



when |cj| < |q| 



(33) 



IV. FERMION DAMPING RATE AT ZERO TEMPERATURE 



In this section, we calculate the fermion damping rate at zero temperature. We first consider the transverse contri- 
bution of gauge field to fermion damping rate. To do this, we will substitute the transverse gauge propagator Eq.((5]) 
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to transvrese self-energy function Eq. tfTTj) . Here, it is convenient to introduce the following spectral representations 

'•+ 00 dun Im[G (wi,k + q)] 



G («e„ + zw m , k + q) = -P / — 

J-00 7T 



1 



|q| 2 + < + n± (iw m ,| q |; 



-P 



+°° dtU 2 1 



— oo 



7T lW m — Ll>2 



-Im 



1 



q| 2 - w I - is s n (^2) 5 + n ± (w 2 , |q|) 



(34) 
(35) 



After carrying out the summation over ui m , we can get 



S T (ie n ,k) = -e 2 ~Tr 



xlm 



1 ' / ttSs* / — Im I G ° ( Wl > k + q )l °j " *%/|q| 2 / — 

J (27T) J-00 7T J-oo 7T 

n s (w 2 ) + n F (uji) 



1 



q| 2 - w 2 - * s S n * + n_L (u) 2 , |q|) 



After analytic continuation ie n — > e + 18, we have 

1 1 



= P- 



1 



ie„ + 0J2 — ui\ s -\- uo'2 — lo\ -f- iS e + L02 — u)\ 
and the imaginary part of fermion self-energy becomes 



— Z7T(5 (e + w 2 — cJi) . 



(36) 



(37) 



ImS T (£,k,T) = -e 5 



rf 2 q 



Im 



|q| 2 - (|k + q| - fi - ef - i<5sgn(|k + q| - fi - e) + Ej.(|k + q| - fi - e, |q[ 



x [n B (|k + q| - (i - e) + n F (\k + q| - /*)] , 
where Eq. (fT?} was used. We now introduce a new variable k' = k + q, and then have 



(38) 



2 f + OO f2Tt 

ImE T (e,k,r) = =■ / d|k'||k'| / dd 

4tt 2 Jo Jo 



xlm 



|k' - k| 2 - (|k'| - n - ef - i(5sgn(|k'| - fj, - e) + n ± (|k'| - /i - e, |k' - k|) 
x [n B (|k'| -n-e) + n F (|k'| - fi)] , 



(39) 



where 9 is the angle between k and k'. Without lose of generality, we suppose that e > 0. 

Now we focus on the zero temperature limit and consider the case at finite temperature in the next section. At 
T = 0, the contribution function reduces to 



n B (|k'| - n - e) + n F (|k'| - ,i) = -6 (|k'| - M ) ( M + £ - |k'|) , 
so the transverse damping rate reduces to 



o2 r» 



Im£ T (e,k) = — / d|k'||k'| / efflm 



2 7; 



1 



|k' - k| 2 - (|k'| - n - e) 2 + iS + II_L(|k'| - n - e, |k' - k|) 



(40) 



(41) 



In general, there are two kinds of approximations: on-shell approximation and fixed-momentum approximation. We 
now consider the on-shell approximation 



e = £ k = £ k - /x = |k| - fi, 



and convert the damping rate to 
2 /-/i+ik 



2 fH+tk f2TT 

ims T (ek) = ^y d|k ' l|k Vo deim 



|k' - k| 2 - (|k'| - |k|) 2 +iS + n ± (|k'| - |k|, |k' - k| 



(42) 



(43) 



The fixed momentum approximation will be discussed later. 



r 



To proceed, we will substitute the analytical expres- last section into this formula. At the T — limit, the 
sion of polarization function Il_i_(u;, |q|) obtained in the 
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integration over parameter x in Eq. ([32|) and Eq. (l33t can 
be analytically carried out. The expression for n_L(w, |q|) 
at T — is presented in the Appendix. Such expression 
is clearly too complicated to be used. In order to get an- 
alytical results for fermion damping rate, it is necessary 
to make proper approximations to TL±(uj, |q|). 

In the present problem, it is important to observe 
that the dominant contribution of the above integral 
comes form the region \w\ <C |q| and |q| <C /i in 
TL±(lu, |q|), so that we can simply the polarization func- 
tions by restricting the energy-momentum to this region. 
This approximation method was used by many authors 
previouslji^££2£. In this region, the polarization func- 
tion can be significantly simplified and is given by 



Ren x (w,|q|) = 



Ne 2 n uj 2 
2tt TqT 



ImllxCw, |q|) ~ -sgn(cj) 



Ne 2 fi H 



(44) 



(45) 



If we take the static limit, uj — > 0, both the real and imag- 
inary parts of the transverse polarization function van- 
ishes, II_l(cj, |q|) — > 0. This implies that the transverse 
gauge field remains massless even after including the dy- 
namical screening effect due to particle-hole excitations. 
This property is robust against higher order corrections 
and indeed a consequence of gauge invariance. However, 
the chemical potential /i does affect the transverse gauge 
interaction between Dirac fermions and thus should af- 
fect the fermion damping rate. Substituting the above 
expressions for H±(uj, |q|) into Eq. (|4"3")l , we finally get 



Im£ T (&)«C(/i)e 



k • 



where 



CV) = -- 



#2. 



(-;?, 



(46) 



(47) 



8\/37r3 N 

Apparently, this damping rate displays non- Fermi liquid 
behavior at zero temperature. 

We now consider the longitudinal contribution to 
fermion damping rate. Starting from Eq.(j4]) and Eq. (|10p 
and then using the same steps presented in the above, we 
write the longitudinal damping rate as 



ImS L (&) 



e 

~4^2 
xlm 



d|k'||k'| 



dO 



|k'-k|2 + n 00 (|k'|-|k 



|k'-k|) 

(48) 

As H±_{uj, |q|), the polarization function IIoo(w, |q|) is 
also too complicated even at T = (see Appendix). In 
the region \u\ <C |q| and |q| <C /i, we have the following 
simplified expressions for IIoo(w, |q|): 

Ne 2 fj, 



Ren 00 (o;, |q|) = 



2tt 



ImIIoo(w, |q|) « sgn(w) 



Ne 2 n M 
2tt TqT' 



(49) 
(50) 



In the static limit, uj 0, the imaginary part 

ImIIoo(a;, |q|) vanishes but the real part Ren o(u;, |q|) is 
a constant. Therefore, from Eq.(j4|), the temporal com- 
ponent of gauge field propagator is found to be 



D m (u = 0,q) 



1 



Ne 2 



(51) 



in the static limit. Comparing with the transverse 
component of gauge field propagator defined by Eq.([S]), 
Eq. (l44"]) and Eq. (|4"5|) . the temporal component has a 
static screening and chemical potential \x defines the De- 
bye screening length. It reflects the effect of particle-hole 
excitations on the initially long-range temporal gauge 
interaction. This is the key difference between Dirac 
fermion systems with zero and finite chemical potential. 
The short-range temporal gauge interaction is expected 
to produce only normal Fermi liquid behavior. Substi- 
tuting them into Eq. (|48|) . it is easy to get 



ImS L (6c) 



(52) 



This expression vanishes faster than £k as £k ~* 
and thus is a normal Fermi liquid behavior. As shown 
ir*2<l, the perturbative result of zero-temperature fermion 
damping rate is divergent at fi — 0. The finite chemical 
potential eliminates the divergence and at the same time 
leads to normal Fermi liquid behavior. 

The total fermion damping rate should be 

ImE(a) = Im£ T (&) + Im£ L (&) « C{^\. (53) 

This result is obtained using the on-shell approximation. 
We can alternatively use the fixed momentum approx- 
imation. The momentum can be chosen as the Fermi 
momentum, so at zero temperature the fermion damp- 
ing rate depends only on the energy e. After explicit 
computation, we found that 



Im£ T (e, |k| = n) 
Im£ L (e, |k| = fx) 



C(/X) £ 3 

1 



2ttA^ 



So the total damping rate is 

Im£ (e,|k| =n) « C(fi)ei. 



(54) 
(55) 

(56) 



This has the same form as that obtained in the on-shell 
approximation, with ^ being replaced by e. 

There are three important features of this damping 
rate. First, when we take the /i — > limit, this result does 
not reduce to the cx e 1 / 2 result obtained at /i = 0. If we 
use the exponent z appearing in the energy-dependence 
e z of damping rate to characterize the ground state of 
the fermion-gauge system, then there is a sudden change 
of ground state once \i departs from zero. It appears that 
the Dirac fermion systems exhibit distinct behaviors at 
zero and finite (i. This difference arises from the differ- 
ence in topology of Fermi surface: at finite [i the system 



9 



has a finite one-dimensional Fermi surface, but at /i = 
the Fermi surface shrinks to a zero-dimensional point. 
Second, at finite fi, as \x grows from certain small value, 
the energy-dependence of fermion damping rate does not 
change. Third, at any fixed energy e the fermion damp- 
ing rate is proportional to /i -1 / 3 , so the Dirac fermions 
become more well-defined as chemical potential grows. 

The fermion damping rate oc e 2 / 3 seems to be a uni- 
versal behavior. It has the same energy-dependence as 
that in two-dimensional non-relativistic fermion-gauge 
systems with a large Fermi surface^—. Such energy- 
dependence also appears in some two-dimensional elec- 
tron systems where fermions interact strongly with fluc- 
tuating ferromagnetic order parameter 24 or fluctuating 
nematic orde r 25 ' 26 , as well as in two-dimensional electron 
systems near a Pomeranchuk instability 27 . 

Using the Kramers-Kronig relation, we get the real 
part of fermion self-energy 



ReS(£) oc v / 3C(/z)sgn(e)e^ 



(57) 



It has the same energy-dependence as the imaginary part. 
It is easy to show that the renormalization factor Z = 0, 
which is the characteristic of a non-Fermi liquid^ 2 .. 

It is also interesting to study the effective QED 3 theory 
without Maxwell term for the gauge field— Now the 
propagator for the gauge field in Matsubara formalism 
has the form 



D 00 (iw m ,q) 



1 



U 00 (iu} m ,q) ' 



Dij (lUJ m ,q) = [Sij - -^y- 



q 2 / II_l (iu m ,q) ' 



(58) 
(59) 



Using this propagator, we found that the longitudinal 
and transverse fermion damping rates are 



ImE T (£k) 
ImS L (^ k ) 



JL 

'ttN 
1 



d 4 



2;rV^ ln (f 



(60) 
(61) 



in the on-shell approximation. In the fixed momentum 
approximation, we have 



ImS T ( £) |k|= M ) « -JL ["dS'^j, (62) 



Im£ L (e, |k| = n) 



^ £ 2 lnf£ 

2TTN/J, \n 



(63) 



Clearly, in both the on-shell and fixed momentum ap- 
proximations, the total fermion damping rate is diver- 
gent. Note that a similar divergence also exists at zero 
chemical potential 2 ^. It seems that such divergences are 
directly related to the absence of Maxwell term for the 
gauge field. 



V. FERMION DAMPING RATE AT FINITE 
TEMPERATURE 

We now consider the fermion damping rate at finite 
temperature. The polarization function at finite T should 
be used when calculating the fermion self-energy. Here it 
will be convenient to adopt an important approximation. 
At low temperature T -C /i, we can still use the polar- 
ization functions obtained at zero temperature. This ap- 
proximation was previously employed in Refs. 3 ' 4 . In the 
limit T <C n, we can simply choose the upper boundary 
value of |k'| as (i + T. The reason is that the fermions are 
primarily scattered into states in the outside of the Fermi 
surface, because most of the states on (and below) the 
Fermi surface are already occupied by other fermions at 
low temperature. The lower limit of |k'| can be assumed 
to be Moreover, at finite temperature, the occupation 
number functions can be well simplified as 



n B (|k'|-|k|)+n F (|k'|-M) 



T 



Ik'|-|k| 



(64) 



After straightforward computation, we finally have 



Im£ T (T) 



12\/37r§ Ni iA 



dS'-r, 

3 



(65) 



which is divergent. It is interesting to note that this 
divergence is very similar to that appearing in the non- 
relativistic fermion-gauge problem (see paper of Lee and 
Nagaosa in Ref^). The longitudinal contribution to 
fermion damping rate at finite temperature is found to 
behave as 



T 2 (T 
ImS L (T) cx — In - 
A 4 \M 



(66) 



which is the typical behavior of normal Fermi liquid in 
two spatial dimensions. Apparently, the total fermion 
damping rate is divergent. 



VI. SUMMARY AND DISCUSSION 

In summary, we studied the effect of finite chemical po- 
tential [i on the damping rate of massless Dirac fermions 
in QED3. At zero temperature, the total damping rate 
behaves as ImX(e,/i) oc /i~ 1 / 3 e 2 / 3 , which vanishes slower 
than e does near the Fermi surface. This non-Fermi liquid 
behavior is primarily generated by the long-range trans- 
verse gauge interaction, while the longitudinal gauge in- 
teraction becomes short-ranged and thus only leads to 
normal Fermi liquid behavior. It is important to note 
that the expression of ImE(e) at fi = can not be ob- 
tained by simply taking the /i — > limit from ImS(e,//) 
at finite [i. This indicates that the fermion damping rate 
displays different e-dependence at zero and finite chemi- 
cal potential, although Fermi liquid theory breaks down 
in both cases. 
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At high fermion density, the Fermi surface becomes 
very large. Now the massless Dirac fermion with linear 
energy spectrum is no longer a good description for the 
low-energy excitations. The system is then described by 
the non-relativistic fermion-gauge theory 4 ^—. Therefore, 
the results obtained in this paper are valid only when fi 
is not too large. 

We have to admit that it is unclear how to get a phys- 
ically meaningful fermion damping rate at finite temper- 
ature and finite chemical potential. When T <C n, al- 
though the longitudinal component of damping rate has 
a normal Fermi liquid result, the transverse component 
ImEx(T) is divergent. At present, there seems to be no 
efficient way to cure such divergence 4 ' 28 . In principle, 
it is possible to get a divergence-free damping rate by 
studying the self-consistent, Eliashberg type, equations of 
fermion self-energy function and polarization functions at 
finite temperature, as we have done previously 20 . How- 
ever, unlike in the case of zero chemical potential—, we 
found it difficult to obtain satisfactory results from the 
corresponding Eliashberg equations at finite chemical po- 
tential. These problems surely deserve more thorough 
investigations in the future. 

Finaly, we also calculated the fermion damping rate 
when the QED3 action has no Maxwell term for the 
gauge field. A divergence appears once the Maxwell term 
is dropped. This divergence has different origin with 
that appearing in the damping rate at finite tempera- 
ture and finite chemical potential, and arises due to the 
absence of Maxwell term. Its appearance may not be 
surprising since we have already met it when studying 
the fermion damping rate at zero chemical potential 20 . 
Unfortunately, we are not aware of any available method 
to eliminate the divergence brought by the absence of 
Maxwell term at both zero and finite chemical potential. 



1. The expression for ImIIoo(w, |q|) 

We first present the expressions for the region |u>| > |q|. 
For < /i < J2tlsL, 



Irnlloo (w, |q|) = sgn(w 



Ne 2 



16 V" 2 -|q| 2 ' 



(Al) 



For 



M-|q| 



^l+|q| 



2 < V < 2 > 
Imlloo (u), |q|) = sgn(w) 



Ae 2 |q| 2 



16tt ^w 2 - |q| 2 



^ A ( 2 m- M 



, (A2) 



where A\{y) = y\/l — y 2 + arcsiny. For \x > M±M 5 

Imn o(^,|q|) =0. (A3) 

We then present the expressions for the region |w| < |q|. 
ForO< M <^ d , 



For JHhH < M < 



Imn 00 (w,|q|) =0. 
|q|+H 



(A4) 



ImUoo (w, |q|) = sgn(o;) 



x.Bi 



Ne 2 |q| 2 



167T y^f^ 
2/i+ \uj\ 



where B x (y) = y^Jy 2 - 1 - In y + yjy 2 - I 



|q|+H 

2 ■ 



(A5) 



For (i > 
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ImII 00 (w, |q|) = sgn(o;) 



5i 



Ne 2 |q| 2 



167T ^q|^ 
2fi+ \uj\ 

T 

2fi - \uj\ 



2. The expression for ReIIoo(w, |q| 



(A6) 



Appendix A: Polarization function at T — 

When calculating the fermion damping rate at zero 
temperature T = in Sec. IV, it is necessary to first 
know the temporal and transverse component of vacuum 
polarization functions. At T = 0, the integration over 
parameter x can be carried out analytically, with the re- 
sults being presented below. Here the chemical potential 
can be taken to be positive: /i > 0. 



We first present the expressions for the region |w| > |q|. 
For < (i < ^dsi, 



Ren 00 (w,|q|) = 



NeV Ae 2 |q| ; 



2n 16tt ^2 _ | q |2 
u\ + 2^ 



Bi 



|q| 

uj\ — 2/i 



(A7) 
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For H_ki < M < M+N ; 



Ren 00 {oj, |q| 

For [i > H±M 

Relloo (w, |q| 



Ne 2 fi Ne 2 



2tt 16tt _ | q |: 



Ae 2 /z Ae 2 |q| ; 



(A8) 



27T 167T y/oj 2 — |q| 2 

1 V |q| 



Si 



2/i- 



|q| 



(A9) 



We then present the expressions for the region \cj\ < |q|. 
ForO<M<^ d , 

Ren 00 (w, |q|) = + 



2ir 16tt v 4qp _ w 2 
2/x- |w| 



(A10) 



For JsbH < M < M+H 



Relloo (u>, |q|) 



Ae 2 /i Ae 2 |q| 2 



+ 



2^ 16tt v 4qp - u 2 

^'(t)]- (a11 > 



For (U > M+H, 



Ren 00 (w, |q[ 



Ae 2 /i 
2tt 



(A12) 



3. The expression for Imllx (a;, |q[) 



We first present the expressions for the region > |q|. 
ForO< M <^ T ^ i , 



Ne 1 

Imn ± (w, |q|) = -sgn(w) — ^co 2 - |q| 2 - (A13) 



For < M < H+N 



Ae 2 



Imn_L (w,|q|) = -sgn(w) J^V^ 2 - |q| 2 

2^ 2 ITT 



(A14) 



where A 2 (y) = —yy/l - y 2 + arcsiny. For ^ > M±lal, 

Imn_L (u, |q|) = 0. (A15) 

We then present the expressions for the region \uj\ < |q|. 
fbr 0< < M^H 



For M=M < M < 



Imn_L (w, |q|) = 0. 
Isl+H 



(A16) 



Ae 2 



Imn± (w, |q|) = -sgn (uj) j^- ^ \q\ 2 - uj 2 



f 2n+ |c 

V |q| 



where B 2 (y) = y^/y 2 - 1 + \n y + \/y 2 -\ 



lal±H 

2 - 



(A17) 



For /x > 



Ae 2 



Imhi (w, |q|) = -sgn (w) Vlql 2- ^ 2 



/2/x+M 

V |q| 

2i< - |w| 



4. The expression for Reflx (lu, |q|) 



(A18) 



We first present the expressions for the region |w| > |q|. 
ForO< A1 <^ T ^, 

Kn , , h Ae 2 M w 2 Ae 2 ^ 

Ren ± ( W ,|q|) = _ —^-M' 



S 2 
B 2 



a;| + 2/A 
\lu\ - 2fi 



(A19) 



For 



M-|q| 



< n < 



H+|q| 



Ae 2 u cj 2 Ae 2 , 



So 



+ 2/A 



|q| / 



(A20) 



For jj, > 



Ae 2 u cj 2 Ae 2 , 



s 2 pM + M 
2 V |q| 

I |q| 



(A21) 
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We then present the expressions for the region \u>\ < |q|. 
ForO< M <^ M , 



(A23) 



ReLU (a,, |q[ 



A^e 2 /x w 2 Ne 



-An 



For < M < M+H 



iVe 2 . 
16^ 


2it + 


m 


|q 









V\<i\ 2 -" 2 



(A22) 



Ren ± (o;,|q|) = + —,/lql 2 



For /i > M+H 



Rellj_ (w, |q| 



ATe 2 ^ w 2 
"^T|qF 



(A24) 
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